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Abstract. The logical calculus for real number system produces new real numbers with 
new operations. New operations naturally produce new equations and thus expand the 
traditional mathematical models which are selected to describe various scientific rules. 
So new operations can describe complex scientific rules which are difficult described by 
traditional equations and have an enormous application potential. As to the equations 
including new operations, engineering computation often need the approximate solutions 
reflecting an intuitive order relation and equivalence relation. However, the order rela- 
tion and equivalence relation of real numbers are not as intuitive as those of decimal 
expansions. Thus, this paper introduces numerical computations to approximate all real 
numbers with decimal expansions. 



1. Introduction 

In 2011, we denned a logical calculus to denote number systems [T]. The number systems 
by logical calculus show typical features as follows: 

1. The logical calculus provides a uniform frame for arithmetic axioms. Based on the 
same logical calculus, small number system can import new axioms to produce big number 
systems. 

2. Consistent binary relation are the nature of number systems. The order relation and 
equivalence relation of each number system is consistent, so all numbers of each number 
system are layed in fixed positions of a number line. 

3. The number systems unify the numbers and operations as a whole, and the only num- 
ber 1 and various operations compose all numbers. So operation distinguishes different 
number systems and is the foremost component of number system. 

4. The logical calculus for real number system produces new real numbers with new 
operations. While producing new real numbers in arithmetic, the new operations certainly 
produce new equations and inequalities in algebra. So the logical calculus not only extends 
real number system, but also extends equations and inequalities. 

In conclusion, the number systems by logical calculus forms a new arithmetic axiom. 
Unless otherwise specified, the number systems, numbers and operations in this note are 
denoted by the logical calculus. Since operations distinguish all number systems, we can 
furtherly define new operations according to the definition of number systems, as is shown 
in ITable II In each row, the column 'Number System' stands for the name of number 
system, the column 'Operation' stands for the name of operations, the column 'Operators' 
stands for the operators of some operation. 

Real operations naturally produce new equations such as y = [x + + + +[1 + 1]], 
y = [[1 + 1] x], y = [x////[l + l]] and so on. In other words, real operations expand 
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Table 1 . Definitions Of New Operations. 



Number System 


Operation 


Operators 


natural number system 


natural operation 


+ 


integral number system 


integral operations 


+,- 


rational number system 


rational operations 


+,++,-,- 


real number system 


real operations 


+,++, + + +, + + ++,••• , 

/J/,///,////,--' 



the traditional mathematical models which are selected to describe various scientific rules. 
So real operations can describe complex scientific rules which are difficult described by 
traditional equations and have an enormous application potential. 

As to the equations including real operations, engineering computation often need the 
approximate solutions reflecting an intuitive order relation and equivalence relation. Al- 
though the order relation and equivalence relation of real numbers are consistent, they are 
not as intuitive as those of decimal expansions. In practice, it is quicker to determine the 
order relation and equivalence relation of decimal expansions. So we introduce numerical 
computations to approximate real numbers with decimal expansions. 

[U Theorem 2.3] implies that decimal expansions are dense in the real numbers with 
the usual topology. So in numerical computations for real number system, all operands 
and outputs are denoted by decimal expansions to intuitively show the order relation and 
equivalence relation. 

The paper is organized as follows. In Section 2, we revise real number system to 
facilitate numerical computations. In Section 3, we introduce the root-finding methods 
for a special equation. In Section 4, we introduce the numerical computations for real 
number system. 

2. Revision For Real Number System 

The following is the revision for real number system, which can facilitate numerical 
computations.. 

Definition 2.1. Real number system is a logical calculus -R{$, ^} such that: 
${ 

(2.1) V{0, a, 6, c, d, e, /, g, h, i, j, k, /}, 

(2.2) C7{0,1, +,[,], -,/,T, _!_,_}, 

(2.3) P{0,G,C,^,|,=,<,<,||}, 

(2.4) V o C{0, a, 6 ■ ■ ■ , 1, + • • • , aa, ab - ■ ■ , al, a + • • • , 6a, 66 • • • , 61, b + • • • , 
aaa, aab ■ ■ ■ , aal, aa + ■ ■ ■ , baa, bab ■ ■ • , 6al, 6a + ■ • • }, 

(2.5) C o C{0, 1, + • ■ ■ , 11, 1 + ■ • • , 111, 11 + • • ■ }, 

(2.6) V o C o P{0, a, 6 • • • , 1, + • • • , G, C • • • , aa, ab - ■ ■ , al, a + • • • , a G, a C • • • , 
6a, 66 • • • , 61, 6 + • • • , 6 G, 6 C • • • , aaa, aa6 • • • , aal, aa + • • • ,aa^,aa<Z--- , 
6aa, 6a6 • • ■ , 6al, 6a + • ■ ■ , 6a G, 6a C ■ ■ • }, 
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(2.16) 
(2.17) 
(2.18) 
(2.19) 



a e V) {{a = 0) V (a = a) V (a = b) ■ ■ ■ ), 
a e C) ((a = 0) V (a = 1) V (a = +) ■ ■ ■ ), 

a G (V o C)) O ((a = 0) V (a = a) V (a = b) ■ ■ ■ V (a = 1) ■ ■ ■ V (a = aa) 
V(a = ab) ■ ■ ■ V (a = al) ■ ■ - V (a = aaa) V (a = aao) ■ ■ • V (a = aal) ■ ■ ■ ), 
(a G (C o C)) ^ ((a = 0) V (a = 1) V (a = +) • • • V (a = 11) V (a = 1+) • • • 
V(a= 111) V (a = 11+)- ••), 

(a G (V o C o P)) ((a = 0) V (a = a) V (a = 6) • ■ ■ V (a =G) ■ ■ • V (a = aa) 
V(a = afe) • • • V (a = a G) • • • V (a = aaa) V (a = aao) ■ ■ ■ \/ (a = aa ■ ■ ■), 
(a G (V o O) A (6 G (V o C)) A(c6(Vo C)) A ( J G (V o C)) A (e G {V o C)) 
A(f e{Vo O) A (g G (V o C)) A (a G (V o C)) A (i G (V o C)) A 

G (y o C)) A (jfc G (V o C)) A (f G (y o C)) A (a G (V o C o P)) A 
1 G (y o C o P)) A (c G (y o C o P)), 
'a C {6, c}) <£> ((a C 6) V (a C c))) A 
a C {&, c, J}) ((a C 6) V (a C c) V (a C J))) A 
a C {&, c, J, e}) <^> ((a C 6) V (a C c) V (a C J) V (a C e))) A 
a C {o, c, J, e, /}) ((a C 6) V (a C c) V (a C d) V 



a C e) V (a C /))) A 
a C {o, c, d, e, f, g}) ((a C 6) V (a C c) V (a C d) V (a C e) V 



V(SCff)))A 

i C {6, c, d, e, /, a}) ((a C 6) V (a C c) V 
V(a C§)V(«C a))) A 
((a C {b, c, d, e, /, g, h, i}) ((a C 6) 
V(a Cj)v(fla)v(aC i))) A 



a C d) V (a C e) V (a C /) 
V (a C c) V (a C d) V (a C e) V (a C /) 



, a 



, a 



C {6, c, d, e, /, g, h, i, j}) ((a C 6) V (a C c) V (a C d) V (a C e) V 
C /) V (a C g) V (a C h) V (a C i) V (a C j))), 
C {b,c,d,e, f,g,h,i,j,k}) ^ ((a C b) 



V (a C c) V (a C d) V (a C e) V 
flC/)V(aCj)V(aa)v(oCi)v(aCj)v(aC fc))), 
a C {o, c, d, e, /, (7, a, ij, k, I}) <^> ((a C b) V (a C c) V (a C d) V (a C e) V 
"C/]v(ay)V(aa)v(ag)v(aCj)v(aa)v(flC [))), 
= dfee/^) A -1(6 C {a, c, d, e, /, (?}) A -.(/ C {a, 6, c, d, e, (?}) A 
->■ a) ||(/ ->■ i)) ^> (aac = dheig), 
abode = f) A -.(6 C {a, c, d, e, /}) A -.(d C {a, 6, c, e, /}) A ((6 ^ ^) || (d -)• a)) 
=^ (a(?cae = /), 
a — )■ 1 1 [a6a] , 
b -»■ +|-, 
c|d -> e|/|a, 



a 
a6c 

(6 



e 



e, 
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(2.20) 
(2.21) 
(2.22) 
(2.23) 
(2.24) 
(2.25) 
(2.26) 
(2.27) 
(2.28) 
(2.29) 
(2.30) 
(2.31) 
(2.32) 
(2.33) 
(2.34) 
(2.35) 
(2.36) 
(2.37) 
(2.38) 
(2.39) 
(2.40) 
(2.41) 
(2.42) 
(2.43) 
(2.44) 
(2.45) 
(2.46) 
(2.47) 
(2.48) 
(2.49) 
(2.50) 
(2.51) 
(2.52) 
(2.53) 
(2.54) 
(2.55) 
(2.56) 



f^-\-f, 

9->/\/g, 
& ->+)!!(<-►-), 

h -t +h)\\(i -> -z), 
i->-)||(ft->+), 

h^+)\\(j^/), 
h^+h)\\(j^/j), 
k-> [l + l]|[l + fc], 

a < [1 + a] , 

a < b) A c (([a + c] < [6 + c]) A ([a - c] < [6 - c]) A ([c - 6] < 
[1 - 1] < a) A (a < 6) A ([1 - 1] < c) =^ ([a - -c] < [6 - -c]), 
[1 - 1] < a) A (a < 6) A ([1 - 1] < c) =>- ([c - -6] < [c - -a]), 



([1-1]<&)=>(1< [a --/&]), 
(K6)^([l-l]<[a/^]), 
(a < b) (1 < [o/ga]), 
(a < 6) A ([1 - 1] < c) =J> ([aec] < [oec]), 
(1 < 6) A ([1 — 1] < c) A ([aec] < [Sec]) (a < b), 
([1 - 1] < 6) A (6 < c) ([aeo] < [aec]), 
([1 - 1] < 6) A ([1 — 1] < c) A ([aeo] < [aec]) =>• (6 < c) 
b] = [a/5]), 

-b] = [a/ mi 



1 < a) A 
1 < a) A 
1 < a) A 
1 < a) A 
1 < a) A 
1 < a) A 
1 < a) A 
a A b =>- ( [a 

SA-i(6= [1- 1]) ([a- 
a =>- ([a - a] = [1 - 1]), 
a A 6 =>- ([a + 6] = [6 + a]) 
a A 6 =>• ([[a - ft] + b] — a) 
a A 6 A c ([[a -b]+c 
a A & A c ([a + [6 + c] 
a A 5 A c =>• ([a + [b - c] 
a A b A c ([a - [b + c] 
a A 6 A c ([a - [6 - c] 
a =^ ([a + +1] = a) A ([a - 
-■(a = [1 - 1]) => ([a a 



s] - 51), 

o] + c]), 
&]-c]), 

a - 6] + c]), 

•i]=a), 
= i), 



a A 6 => ([a + +6] = [b + +a]), 
aA6Ac^([o + +[6 + +c]] = [[a + +6] + +c]), 
a A 6 A c ([a + +[6 + c]] = [[a + +6] + [a + +c]]), 
a A 6 A c ([a + +[& - c]] = [[a + +6] - [a + +c]]), 
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a A 


-,(6 = fl - 11) (ffa - -61 + +61 = a), 

\ L J / V LL J J / * 






(2.58) 


a A 


-1(6 = fl — 11) A c =>- ((ffa 61 + +c1 = ffa + +c1 - 

V L J / VVLL J 1 1 J LL 1 1 J 


-61) A 

J / 






(ffa 

V LL 


+ c] - -61 = ffa - -61 + fc - -611) A (ffa — cl — -61 

1 J J LL J 1 L JJ/ VLL J J 








ffa- 

LL 


- -61 - fc - -611)), 

J L J J / / * 






(2.59) 


a A 


-.(6 = fl - 11) A -.(c = fl - 11) ((fa + + f6 - -cl] = 

\ L J/ V L J / \ V L 1 1 L JJ 


= ffa + +61 - 

LL 1 1 J 


-cl) A 

J / 




(fa- 


- -f6 + +c]] = ffa - -61 - -cl) A ([a - -\b - -cl] = 

L 1 1 J J LL J J/VL L JJ 


ffa - -6] + 

L L J 1 


+cl)), 

1 J / / * 


(2.60) 


a =r> 


(la + + + 11 = a) A (la 11= a), 

V L 111 j / V L J / ? 






(2.61) 


a => 


(\l + + + a] = 1), 

V L 1 1 1 J / ' 






(2.62) 


-•(a 


= fl - 11) (fa + + + fl - 111 = 1), 

L J / V L 1 1 1 L JJ / J 






(2.63) 


([1 - 


- l] < a) (ffl - 1] + + + a] = [1 - 11), 

J / VLL J III J L J/" 






(2 64) 


(fl - 


- 11 < a) A -1(6 = fl - 111 A c ((11a 61 + + 4 


- 61 = a) A 






([[« 


6] + + + cl = [fa + + + cl 61) A (fa + + 

J JLL J J/\L 


+ [c - -6]] 


= 




f[a- 


h + + c] 6])), 






(2.65) 


([1- 


- 1] < a) A (fl - 1] < 6) A c (([a + + + [b///a\] = 


6) A 






([[« 


+ + + c]///b] = [c + +[a///b]}) A ([[a - -6] + + + c 


] = 






[fa- 


h + + c] - -[6 + + + c]])), 






(2.66) 


([1- 


- 1] < a) A (fl - 1] < 6) A (fl - 1] < c) =► (([[a///c] 


--iv/m 





(2.67) 



A ([[a + +6]///c] = f[a///c] + [6///c]]) A ([[a - -6]///c] = 
- [&///S]])), 

.(a = [1 - 1]) A -.(6 = [1 - 1]) A -.(c = [1 - 1]) ^ (([[a + +6] + + + c] = 
a + + + c] + +[6 + + + c}}) A (fa + + + [6 + +c]] = [fa + + + 6] + + + c]) 
A([a + + + [6 + c]] = [fa + + + 6] + +[a + + + c]]) A ( [a + + + [6 — c]] = 





f[a + 


+ + 6] [a + + + c]])), 






(2.68) 


(fl- 


1] <a)A-(6= [l-l])A-.(c= [1- 


-1]) 


=►(([5 [6 + +c]] 




[fa- 


- -6] c]) A (fa [6- - 


c]] = 


[fa 6] + + + c])), 


(2.69) 


([1- 


1] < a) A (fl - 1] < 6) A^(c= [1 - 


1]) = 


K([[a + +6]- -c] = 




[fa — 


--c] + +[6 c}}) A ([[a --6] 




-s] = 




[fa — 


~-c]--[b c]])), 






(2.70) 


(1 < 


a) (fa + el] = a), 






(2.71) 


(1 < 


a) (fa + +e[l- 1]] = 1), 






(2.72) 


(1 < 


a) (fa - /l] = a), 






(2.73) 


([1- 


1] < a) (fl + +ea] = 1), 






(2.74) 


([1- 


1] < a) (fl - -fa] = 1), 






(2.75) 


(! < 


a) =► ([l//ga] = [1 - 1]), 






(2.76) 


(! < 


a) = 1), 






(2.77) 


(1 < 


a) A (fl - 1] < 6) ([[aib]hb] = a), 






(2.78) 


(1 < 


a) A (fl - 1] < 6) =>■ ([[ahb]ib] = a), 
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[2.79) 


(1 < a) A (1 < 6) 


(\bh\ajb}} = 


a), 




[2.80) 


(1 < a) A (1 < 6) 


([[6M7&1 = 

VLL \J \ 


a), 




'2.81) 


(1< a) A (1< b) 


(fa + e&l = 

V L 1 J 


[ae[a + e[o - llll) , 

L L 1 L J J J / ' 




[2.82) 


T_l_l_= [1-1], 








[2.83) 


T 1 [1 + 11 =1, 

— — L 1 J — 1 








[2.84) 


^_1_1_= 1, 








[2.85) 


_L 1 [1 + 11 =1 








[2.86) 


T k + 11 + +1] 

— — LLL 1 J 1 1 J 


-11 = T 


[A; -11 / , 

-L J — — ? 




[2.87) 


-LJUK 1 + 1] + + l ] 


- 1]_ = ±_ 


.[fc - l]J_, 




[2.88) 


T_A;_[[1 + 1] + +/]_ 


- = ru* - 


l]J_ + T_[fc-l]_ 


_P + 1]J, 


[2.89) 


-LJU[1 + 1] + 


- = - 


1]J_+-L_[A;-1]_ 




[2.90) 


(1 < a) A (T_6_c_) 


A (l__6_c_ 


) ([a + /i[T_6_c_ 


_ - 




[[a + hT_b_c_] - L 


L_5_cJ) 







}, 

2.91) 

2.92) (a C 6) <£> (6 = cad), 

2.93) (a -» feed) A (c -)• e) =>- (a ->■ bed), 

2.94) -»■ o|c) =^ ((a ->• 6) A (a -»■ c)), 

2.95) (a|5 ->■ c) ((a -> c) A (6 -> c)), 

2.96) (a < 6) -1(6 < a), 

2.97) (a < 6) -.(a = 6), 

2.98) (a < 6) A (6 < c) =>• (a < c), 

2.99) (a < 6) A (a G (C o C)) A (6 G (C o C)) 4(aA 6), 

2.100) (a < feed) A (c = e) (a < bed), 

2.101) (afc < d) A (6 = e) =>- (aec < d), 

2.102) (a < feed) A (c — >■ e) A -.(c C {a,o,d}) (a < bed), 

2.103) (a < fecdee) A (c -> /) A -.(c C {a, 6, d, e}) (a < fr/d/e), 

2.104) (a6c < d) A (b -> e) A -1(6 C {a, c, J}) =>- (aec < J), 

2.105) (aac < dbe) A (6 ->■ f) A ->(& C {a, c, J, e}) =>• (a/c < dfe), 

2.106) (aac < dbebf) A (0 — > g) A ->(& C {a, c, d, e, /}) (age < dgegf), 

2.107) (aacW < e) A (6 ->■ /) A -1(6 C {a, c, d, e}) =>• (afcfd < e), 

2.108) (abebd < ebf) A (0 — >■ #) A ->(6 C {a, c, d, e, /}) =>- (agegd < eg]), 

2.109) (aacad < A (6 ->■ a) A ->(& C {a, c, d, e, /, (?}) ^> (ahchd < ehfhg), 

2.110) a = a, 

2.111) (a = 6) (6 = a), 

2.112) (a = 6) -.(a < 6), 
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(2.113) (a = bed) A (c = e) =>• (a = feed), 

(2.114) (afec) A (6 = d) (afec = adc), 

(2.115) (a = bed) A (c — >■ e) A ->(c C {a, fe, d}) (a = feed), 

(2.116) (a = fecdee) A (c — >• /) A -i(c C {a, fe, d, e}) =>■ (a = fe/d/e), 

(2.117) (afec = dfee) A (fe -»■ /) A ->(& C {a, c, d, e}) (a/c = d/e), 

(2.118) (afec = dbebf) A (6 -)■ g) A ->(& C {a, c, d, e, /}) =>- (age = dgegf), 

(2.119) (afecfed = ebfbg) A (fe ->■ a) A -i(fe C {a, c, d, e, /, #}) =S> (ahchd = ehfhg), 

(2.120) (a < fe) ^ ((a < fe) V (a = fe)), 

(2.121) (a < fe) A (fe < c) (a < c) 

}• 

The axiom [TJ (4.125)] is expanded to be the axiom (12. 3p . The axioms [TJ (4.153)]~[TJ 
(4.155)] are expanded to be the axioms (j2.3ip ~ (12.341) . The axiom f!2.35j) is added to the 
real number system. The axioms [TJ (4.157)]~[TJ (4.160)] are expanded to be the axioms 
([2T37j) ~ fl2^0l . The axiom pj (4.161)] is deleted. The axiom pfl (4.167)] is revised to be 
the axiom (I2.46p . The axiom [TJ (4.172)] is expanded to be the axiom (12.511) . The axiom 
P, (4.179)] is revised to be the axiom (12.581) . The axiom [TJ (4.181)] is expanded to be 
the axiom (12.601) . The axiom [TJ (4.185)] is revised to be the axiom f)2.64p . The axioms 
[TJ (4.191)]~P (4.195)] are expanded to be the axioms (gZDD ~ (gZZ5D - The axioms [JJ 
(4.197)]~[TJ (4.200)] are expanded to be the axioms fl2TTTp ~ fl2~m The axiom [JJ (4.201)] 
is expanded to be the axiom (12.811) . The axiom [JJ, (4.210)] is expanded to be the axiom 
(I2.90p . The axioms (l2.120pMI2.12ip are added to the real number system. 

It should be noted that fl2~82l ~ f[2T90l) restrict [T_fe_c_ -L_fe_c_] to be a Farey 

fraction. In the following, we will deduce some numbers and equalities as examples. 



(Al) 


(a -> 1 \[aba\) =^ (a ->• 1) 


fewfl2.16p.d2.94p 


(A2) 


=>- (a — >■ [aba]) 


few (12. 16p. 02.941) 


(A3) 


=>- (a ->• [[afea]fea]) 


by(A2) : Q2.93P 


(AA) 


=► (a [[1 - 1] - 1]) 


6y(Al).(|2.17|) 


(A5) 


(a < [1 + a]) (1< [1 + 1]) 


fewrt2.30p.fAlUI2.105p 


(AQ) 


1 


few()2.99|) 


(A7) 


=► [1 + 1] 




(A8) 


([aba] < [1 + [aba]]) 


fey(A5),(A2),(|2JH5D 


(A9) 


=► ([1 + 1] < [1 + [1 + 1]]) 


6wfAlTfl2.17p.fl2.105p 


(A10) 


=► ([1 - 1] < [1 + [1 - 1]]) 


bv(A8). (Al). (12.171). (12.1051) 


(All) 


=► ([1 - 1] < [[1 - 1] + 1]) 


fewfl2.44p.fl2.100p 


(A12) 


=► ([1 - 1] < 1) 


few(l2.45p.fl2.100P 


(A13) 


=► ([1 - 1] < [1 + [1 + 1]]) 


by(Al2),(A5),(A9),<^m 


(AU) 


=► [1 - 1] 


fey(Ai2),(j2.99|) 


(Alb) 


=>[1 + [1 + 1]] 


by(A13),UM 


(A1Q) 


=^ ([1-1] < [1 + 1]) 


by(A12),(A5) J UM> 
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(All) =*([[1 + 1]- -[1 + 1]]<[[1 + [1 + 1]]- -[1 + 1]]) by(A16),(A9),(^m 

(Ai8) =>- (i< [[i + [i + 1]] - -[i + 1]]) by$n^,$nm$ 

(419) =*[[1 + [1 + 1]] --[1 + 1]] byUM 

Then we deduce the numbers from 

[1 - 1], [[1 - 1] - [[1 + 1] [1 + 1]]], [[1 + [1 + 1]] [1 + 1]] • • • 

(51) (T_k_[[l + 1] + +/]_ = [T_[k - l]_l_ + T_[A; - l]_[l + 1] J) by$ZM) 

(52) (T_[l + 1]_[[1 + 1] + +/]_ = [T_[[l + 1] - 

T_[[l + 1] - l]_[l + 1]J) MM, (EUD, 

(53) =► (T_[l + 1]_[[1 + 1] + +1]_ = [T_[[l + 1] - 1]_1_+ 

T_[[l + 1] - 1]_[1 + 1] J) by^M, (M), 

m =► (t_[i + i]_[[i + 1] + +i]_ = [t_[[i - 1] + i]_i_+ 

t_[[i - 1] + i]_[i + 1] J) TO 

(55) (T_[l + 1]_[[1 + 1] + +1]_ = [T_l_l_ + T_l_[l + 1] J) 6y(|235), ( jZHBj) 

(56) (T_[l + 1]_[1 + 1]_ = [T_1_1_+T_1_[1 + 1]J) 6v(|23Ijl. (ErTTIl i. 

(12TT3D 

(57) => (T_[l + 1]_[1 + 1]_ = [[1 - 1] + 1]) 6y(I2S2I), (EH3J, 

f l2TT3|) 

(58) =>> (T_[l + 1]_[1 + 1]_ = 1) 6^(123511,(12113} 

(59) ([1 - 1] < T_[l + 1]_[1 + 1]_) by{A12),$nM 
(510) (T_[l + 1]_[1 + 1] J byQZMb 
(fill) (±_k_[[l + 1] + +/]_= [-L_[fc-l]J_+±_[fc-l]_[Z + l]J) 6y(EH9D 

(512) (i__[l + 1]_[[1 + 1] + +l}_ = [i__[[l + 1] - 1]_Z_+ 

-!_[[! + 1] - + i]J) 6yd22ED, (EH, 

( 12+181) 

(513) =► (±_[1 + 1]_[[1 + 1] + +1]_ = [±_[[1 + 1] - 1]_1_+ 

+ 1] - i]_[i + i]j) M223}> ra>» 

( 12+181) 

(514) =► (+_[1 + 1]_[[1 + 1] + +1]_ = [±_[[1 - 1] + 1]_1_+ 

-L_[[l - 1] + + 1]J) by®Z®, (E+H 

(515) (±_[1 + 1]_[[1 + 1] + +1]_ = [_L_1_1_ + +_1_[1 + 1]J) 6y(l235D, ( 12+131) 

(516) (+_[1 + 1]_[1 + 1]_ = [+_1_1_++_1_[1 + 1]_]) 6y(jmD,(j2+n]), 

(123133 

(517) =4> (J__[l + !]_[! + 1]_= [1 + 1]) by^M,UM 
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(518) ([1 - 1] < ±_[1 + 1]_[1 + 1]J by(A16),$ZM% 

(519) (±_[1 + 1]_[1 + 1] J by^Mj 

(520) ([[l + l] + e[[l + [l + l]]--[l + l]]] = 

[[1 + l]e[[l + 1] + e[[[l + [1 + 1]] - -[1 + 1]] - 1]]]) by(A5) : (A18), 

(Ai6),d2IEU 

(521) =► ([[1 + 1] + +e[[l + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + e[[l + 1] + +e[[[l + [1 + 1]] - -[1 + 1]] - 1]]]) by<^m>, pTTg]) 

(522) =► ([[1 + 1] + + + e[[l + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + +e[[l + 1] + + + e[[[l + [1 + 1]] - -[1 + 1]] - 1]]]) by®W, (EIBJ 

(523) =► ([[1 + 1] + + + +[[1 + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + + + [[1 + 1] + + + +[[[1 + [1 + 1]] - -[1 + 1]] - 1]]]) by^M, (EHBJ 

(524) =► ([[1 + 1] + + + +[[1 + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + + + [[1 + 1] + + + + 

[[[i - -[i + i]] + [[i + i] - -[i + 1]]] - i]]]) (Em 

(525) =► ([[1 + 1] + + + +[[1 + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + + + [[1 + 1] + + + +[[[1 - -[1 + 1]] + 1] - 1]]]) 6y(A16), d232]), 

(526) =► ([[1 + 1] + + + +[[1 + [1 + 1]] - -[1 + 1]]] = 

[[i + 1] + + + [[i + 1] + + + +[[[i - -[i + 1]] - 1] + 1]]]) by&m, ( PmsD 

(527) =► ([[1 + 1] + + + +[[1 + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + + + [[1 + 1] + + + +[1 - -[1 + 1]]]]) ( I2~TT3|) 

(528) ([[1 + 1] + h[T_[l + 1]_[1 + 1]_ - -+_[1 + 1]_[1 + 1] J] = 

[[[1 + 1] + hT_[l + 1]_[1 + 1]_] - i±_[l + 1]_[1 + 1] J) by(A5), (510), 

(519), ([210]) 

(529) ([[1 + 1] + h[l - -[1 + 1]]] = [[[1 + 1] + hi] - i[l + 1]]) M58), (517), 

(EZHD 

(530) =>> ([[1 + 1] + +h[l - -[1 + 1]]] = 

[[[i + 1] + +hi] - + 1]]) fcytra, mM 

(531) ([[1 + 1] + + + /i[l [1 + 1]]] = 

[[[1 + 1] + + + hi] i[l + l]]) 6y(j2+4D,0l23D 

(532) =►([[! + i] + + + + [i __[i + i]]] = 

[[[1 + 1] + + + +1] [1 + 1]]) (E22D 

(533) ([[1 + 1] + el] = [1 + 1]) by(A5),(^m> 

(534) ([[1 + 1] + +el] = [1 + 1]) by(^J^,<mm, 

f l2TT5|) 

(535) =>> ([[1 + 1]+ + + el] = [1 + 1]) ^(J2+9]),(|2+II]), 

f l2+T5D 
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(£36) ([[1 + 1] + + + +1] = [1 + 1]) fos/d2HSI>, d2HTTJ>, 

(EHUD 

(£37) =► ([[1 + 1] + + + +[1 - -[1 + 1]]] = [[1 + 1]- -[1 + 1]]) fy(£32),(£36), 

f l2TT3|) 

(538) =► ([[1 + 1] + + + +[[1 + [1 + 1]] - -[1 + 1]]] = 

[[1 + 1] + + + [[1 + 1] [1 + 1]]]) by(B27), (£37), 

f l2TT3|) 



Then we deduce the equalities on deducible numbers from £{$, ^/}: 

[[1 + 1] + +[[1 + 1] [1 + 1]]] = [[[1 + 1] [1 + 1]] + +[1 + 1]], 

[[1 + !] + + + +[[! + [1 + 1]] --[1 + 1]]] = [[1 + 1] + + + [[1 + 1] [1 + 1]]] 



The deducible numbers correspond to the real numbers as follows: 
[[1 - !]-[[!+ [1 + 1]] [1 + 1]]] = 

* 1 

[[1 - 1] - [[1 + 1] [1 + 1]]] = -#2, 

[[1 - 1] - [[1 + 1] [1 + [1 + 1]]]] = -^2, 

[[1 - 1] - 1] = -1, 



[[1-1] -[!--[! + !]]] = -\ 



[1-1] = 0, 



[i - -[l + 1]] = ~ 



[[1 + !]///[! + [1 + 1]]] = log 3 2, 



1 = 1, 
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[[1 + 1] [1 + [1 + 1]]] = ^ 

[[1 + 1] [1 + 1]] = #2, 

[1 + [1 " "[I + 1]]] = \, 

[[1 + [1 + !]]///[! + 1]] = log 2 3, 

[[1 + [1 + 1]] [1 + 1]] = 

[1 + 1] = 2, 

[[1 + 1] + [1--[1 + 1]]] = | 

[1 + [1 + 1]] = 3, 
[[[1 + [1 + 1]] [1 + 1]] + +[1 + 1]] = 

3. Root-finding Methods For A Special Equation 

Common root-finding methods and their convergence conditions are listed in jU TABLE 
PT2.3]. In this section, we only consider the special equation as follows: 

(1) The equation are continuous on a closed interval [a,b]; 

(2) The equation have only one root on [a, b]; 

When a acts as the lower guess and b acts as the upper guess, both the bisection 
method[U §5.2] and Brent's methodjU §6.4] always converge and find the only root on 
[a,b]. But Brent's method converges faster than the bisection method and thus becomes 
the main root-finding method for above equation. 

4. Numerical Computations For Real Number System 
4.1. Computation Precedence. 

Definition 4.1. In real number system, unary numerical computation takes the number 
1 as the only operand and outputs the approximated decimal expansion 1.0. 



12 



PITH XIE 



Definition 4.2. In real number system, binary numerical computation takes two decimal 
expansions as operands and outputs an approximated decimal expansion. 

According to §2} each real number except 1 is iteratively deduced between '[' and ']' 
from only one real operation and the other two real numbers. So after iterative unary 
numerical computation on given real number, iterative binary numerical computation on 
the given real number can gradually approximate all real numbers between the outmost 
'[' and ']' with decimal expansions. 

In summary, the computation priority is sorted in descending order as follows: 

(1) Unary numerical computations have precedence over binary numerical computa- 
tions. 

(2) Left unary numerical computations have precedence over right unary numerical 
computations. 

(3) Binary numerical computations on two decimal expansions have precedence over 
those involving a non-decimal expansion. 

(4) Left binary numerical computations have precedence over right binary numerical 
computations. 

Unary numerical computation is easy, so we focus on binary numerical computation in 
the following subsections. 

4.2. Division Of Real Operations. According to the complexity of numerical com- 
putations, we divide real operations into low operations, middle operations and high 
operations. ITable 21 lists their elements in detail. 



Table 2. Division Of Real Operations. 





Low Operations 


Middle Operations 


High Operations 


Operators 


+,++,-, ,/,// 


+ + +, ,/// 


+ + ++, + + + + +,-■• , 
///A/////,"--' 



4.3. Binary Numerical Computations For Low Operations. In real number sys- 
tem, / is equal to — while //is equal to . From a traditional viewpoint, the low oper- 
ations "+, ++, — , " are equal to basic arithmetic operations "+, x, — , So binary 

numerical computations for low operations have been achieved in elementary arithmetic. 

4.4. Binary Numerical Computations For Middle Operations. From a traditional 

viewpoint, + + + is an exponentiation operation, is a root-extraction operation 

and /// is a logarithm operation. In this subsection, we import the binary numerical 
computations for the middle operations "+ + +, , // /" in [21 §23]. 

Let e be Euler's number. Suppose that a G (— oo, +oo) is a decimal expansion, n G Z 
and k G N '. The binary numerical computation for [e + + + a] can be achieved with the 
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Taylor-series expansion as follows. 

lim 2^[[a 



e + 



n 



k— v+oo 



[ft!]] 



n=0 



« [1 + [a - -[1!]] + [[a + + + 2] - -[2!]] + [[a + + + 3] - -[3!]] + ■ ■ ■ + 
[[a + + + n] - -{n\}} + •• •] 

Suppose that a G (0, +oo) is a decimal expansion and h £ Z. Let 6 = [[a— 1] [a+1]], 

then the binary numerical computation for [a///e] can be achieved with the Taylor-series 
expansion as follows. 



[a///e] 



[[[l + &]--[l-6]]///e] 







2 + + 


lim I 




k— >+oo I 



( [[b + + + [2n + 1]] - -[2n + l]] 

n.=0 

[2 + +[6 + [[6 + + + 3] - -3] + [[6 + + + 5] - -5] + • ■ ■ + 
[[6 + + + [2n + 1]] --[2n + !]] + ■■■]] 



Suppose that a and b are two decimal expansions, where a G (0, +oo) and b G (— oo, +oo). 
Then the binary numerical computation for [a + + + b] can be divided and conquered with 
the identity [a + + + b] =Je + + + [b + +[a///e]]]. 

Suppose that [a + + + b] is a real number, where a G (— oo, 0] and b G (-co, +oo) are 
two decimal expansions. Then the binary numerical computation for [a + + + b] can always 
be equated with the basic binary numerical computations as above and basic arithmetic 
operations with the axioms (j2.ipMI2.69p . 



4.5. Binary Numerical Computations For High Operations. The floor function 
[•J and mod operator are imported from [31 §1]. Suppose that the constants n, a±, a^, 
bi, &2 ; Ci, C2 G N. Suppose that the constant d G R with [1 — 1] < d. Suppose that the 
constant e G R with 1 < e. 

Zero can be equated with a fraction [[1 — 1] 1]. Any non-zero decimal expansion 

can be equated with a fraction [[1 — 1] ± [a\ a 2 ]] for a\, d 2 G N. 

The Euclidean Algorithm(EA)[31 §4] can compute the greatest common divisor (di, 62). 

If di and d 2 are divided by (61,02), then the fraction [d\ 62] is fully reduced to 

irreducible fraction. So EA suffices to equate any non-zero decimal expansion with an 
irreducible fraction. The irreducible fraction of a fraction d is denoted as 1(d). 

Theorem 4.3. For a function / : R \ (—00, 1) — > R defined by f(x) = [x + + + ed] , it is 

injective continuous and maps the domain [l,+oo) to the range [l,+oo). 

Proof. According to f)2.19p . the symbol 'e' represents some successive '+' — "+•■■+". 
According to (I2.20p . the symbol '/' represents some successive ' — ' — "— According 
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to f)2.2ip . the symbol 'g' represents some successive '/' — "/ ••■/"■ 





titiogp fna^" -i o^n 1 — I— r~v^l ^xnfn ^h-i =rz ^f'n 
Oll|J|JvJoc nidi Uj\) Uj2 ^ IX- | ?k-J J Willi Jb\ -t— Ju^. 

vviijiiCLiL icoo ui J^dici dii ij^y ; wc dooliiiiv; nidi Jb\ \ 




(A?) 

(Sit) 






/ /l Q\ 

[A3) 


[1 — 1J < a 


by (Premise) 


(A4) 


[ii + + + erf] < [x 2 + + + erf] 


by (A2),(A3),(j2.37!) 


[AS) 


/(xi) = [xi + + + erf] 


by (Premise) 


(A6) 


=> /(ii) < [x 2 + + + erf] 


by (A4),(A5),d210lD 


(A7) 


/(5 2 ) = [s 2 + + + erf] 


by (Premise) 


(AS) 


< f{x 2 ) 


by (A6),(A7),(|2100|) 


(A9) 


=► -.(/(ii) = /(i 2 )) 


bv (12.971) 



(A1)~(A9) derive that /(x) is injective. (A1)~(A8) derive that f(x) is strictly increas- 
ing. 

For any number xq G [1, +oo) and any number [1 — 1] < e, we can always construct So 
as follows: 



(51) [l + + + erf] < [x + + + erf] bv ([2T3TD . fl2TT0|) 

(52) =4> 1 < [x + + + erf] bv (127751) . (I21UI]) 

(53) [1 - 1] < [[x + + + erf] - 1] by <^M) 

(54) [1 - 1] < e by (Premise) 

(55) [[1 - 1] - -[1 + 1]] < [e - -[1 + 1]] by (|Z32D 

(B6 ) ^ [[i _ _[i + 1]] _ [i _ _[i + 1]]] < [e _ _[i + 1]] by (E3hd, mm 

(B7) [1-1] < [e--[l + l]] by dUSD, flSTMP 

(58) =*[[i_i] + [ e - -[i + i]]]<[[£--[i + i]] + [£- -[i + i]]] byjzau) 

(59) [[1 - 1] + [£ - -[1 + 1]]] <[[£ + £]- -[1 + 1]] by <gm>, (I2TTU0J) 

(510) [[£ - -[i + 1]] + [i - 1]] <[[£+£]- -[i + 1]] by (hbd, dnnn) 

(511) =►[[[£- -[i + i]] + i]-i] <[[£ + £]- -[i + i]] by (E3HD, (EMD 

(512) [[[£- -[1 + 1]]- 1] + 1] <[[£ + £] --[1 + 1]] by (J23BD, dHOU) 

(513) =^ [£--[1 + 1]] <[[£ + £] --[1 + 1]] by flUSD, d37Tmj) 

(514) =► [e - -[1 + 1]] < [[£ + +[1 + 1]] - -[1 + 1]] by <$M), 

(515) [£--[1 + 1]] < [£ + +[[1 + 1] - -[1 + 1]]] by fl239D,fl2l00D 

(516) =^[£--[l + l]]<[£ + +l] by (&M,(gim 

(517) ^ [£--[1 + 1]] <£ by 02311), J3IDDJ 

(518) S = [£--[1 + 1]] 



We construct S according to [[xo + + + erf] — So]. 
(1) 1< [[xo + + + ed]-8 ]. 
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We construct S as follows: 

(CI) S 1 = [x -[[[x + + + ed)-S ] fd]] 

(C2) S 2 = [[[[x + + + ed]+S ] fd]-x ] 

(C3) S = min{S 1 ,S 2 } 



(Dl) 




r r r-i r -i n 

8 = [e [1 + 1]] 




by (B18) 


{D2) 




=>- [1 - 1J < d 




by 










(Iz.lUUp 


[L>6) 




F i 1 + -"-J J < e 




D Y l-D- 1 ' ) 


(DA) 




=>- <5 < £ 














<\2 l m h 


(D5) 




S = min <5 2 } 




bv (C3) 


(D6) 




S < 5i 






(D7) 










{-DO) 




<?i = [5 - [po + + + erf] - 5 ] /dj] 






( nql 


[Xq 


- 5x] = [x - [x - [[[x + + + ed] - S ] 


fd)]] 


hv 119 1 1 Ah 
uv iii. ii^ji 




Fo 


- S^ = [[x - x ] + [[[x + + + ed] - S ] 


- M] 


DV IIZ.0U|I.IIZ.llO|l 


( ^)^^\ 
[Mil.) 


Fo 


- 5i] = [[[[x + + + ed] -S ] fd] + [x - 


~ x o]] 


L.Tr (19 A A h (19 1 1 'ih 
DV IIZ.44II.IIZ.110II 


( n^ c )\ 


=> [x 


- St] = [[[[[x + + + ed] -So] fd] + x ] 


- Xo] 


Kir (19 /I Ch (19 1 1 
DV (IZ.4oH.IIZ.llo!) 


( m "X\ 

[JJi6) 


Fo 


- S^ = [[[p + + + ed] - So] fd] - x ] + x ] 


DV (IZ.4DP.IIZ.110II 


( m a\ 




=> [x - Sx] = [[[xo + + + ed] - So] fd] 




Lnr (19 A Kt\ (19 1 1 'ih 
DV (IZ.40p.IIZ. 11011 


( m ^ 

{ULO) 




l<[[[xo + + + ed]-S ]- -fd] 




Dy (^rremisejjijZiO^, 










rt9 74h 

IIZ. ( 4]) 






=>- 1 < [x o - <5i] 




Dy ( i U14J,( i JJl0j, 










<\'> 1 0Oh 

IIZ. 1UUII 






=> [x - S^ < [x - S] 




Lar (Tlfi"! (19 'il h 


(018) 




[[xo - Si] + + + ed] < [[x - S] + + + ed] 




bv fD17).|2.37|, 










(12.114[) 


(D19) 




> [[po + + + ed] - So] fd] + + + ed] < 










po ~ S] + + + ed] 




by (Dl4),fl2jrjlj) 


(D20) 




=> po + + + ed] - S ] < p - S] + + + ed] 




bv 02.77p.fl2.101D 


(D21) 




=> po + + + ed] - e] < [[x + + + ed] - S ] 




by (D4),f|2.3ip 


(D22) 




[[xo + + + ed] - e] < [[x - S] + + + ed] 




by (D21),(D20), 










(|2.98|) 


(D23) 




^[f( Xo )-e]<f([x -S]) 




by (Premise) 


(D2A) 




S 2 = [[[[x + + + ed] + S ] fd] - xo] 




by (C2) 
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(D25) 


[s&o + 6 2 ] = [x + [[[[x + + + ed] + S ] fd] - 


x ]] 


by (j2.114|) 


(D26) 


=> [x + S 2 ] = [[[[[x + + + ed] + S ] fd] - x ] 4 


-x ] 


bv fl2.44p.fl2.113p 


(D27) 


[x + 5 2 ] = [po + + + ed] + 60]- - fd] 




bv fl2.45p. (12.1131) 


(D28) 


[1 - 1] < 6 




by (D2) 


(D29) 


[[1 - 1] + 1] < [So + 1] 




by (D28),(j2 i 31]) 


(D30) 


K [So + 1] 




bv (I2.45p.fl2.l0ip 


[JJOL) 


1 <. [1 -|- OoJ 




UV lli.^^jl . IIZ. ±UUII 


(D32) 


1 < [xo + + 4- ed] 




by (B2) 


(L>33) 


[1 + So] < [[$0 + + + ed] + 5o] 




by (]2.31|) 


(L>34) 


1 < [[x + + + ed] + S Q ] 




by (D31),(D33), 
(|2.9«P 


(D35) 


=> 1 < [p + + + ed] +5 ] fd] 




by (]2.34|) 


(D36) 


po + 5 2 ] + + + ed] = 








[[po + + + ed] +5 ] fd] + + + ed] 




by (D27),(D35), 








(]2.114| 


(D37) 


=> Po + S 2 ] + + + ed] = p + + + ed] + S ] 




bv fl2.77p.fl2.113D 


(D38) 


5 <e 




by (D4) 


(D39) 


[S + [xq + + + erf]] < [e + [£q + + + ed]] 




by (|2.3ip 


(L>40) 


=> po + + + ed] + 5 ] < [s + [x + + + ed]] 




bv O2.44p.fl2.10ip 


(D41) 


=> p + + + ed] + So] < p + + + ed] + e] 




bv fl2.44p.fl2.100p 


(D42) 


=^ po + S 2 ] + + + ed] < p + + + ed] + e] 




by (D41),(D37), 








(]2.10l|) 


(D43) 


po + + + ed] - 5 ] < [[x + + + ed] - [1 - 1]] 




by (D2),(P3U 


(L>44) 


=> po + + + ed] - So] < [[[x + + + ed] - 1] + 1] 




bv fl2.50p.fl2.100D 


(D45) 


po + + + ed] - <J ] < [i + + + ed] 




bv fl2.45p.fl2.100p 


(D46) 


po + + + ed] - S ] = 








[[po + + + ed] - So] fd] + + + ed] 




by (Premise), (J2.77P 


(D47) 


=> [x + + + ed] = 








[po + + + ed] /d] + + + ed] 




by (B2),(J2.77P 


(D48) 


[[Po + + + ed] - 5 ] fd) + + + ed] < 








[po + + + ed] fd] + + + ed] 




by (D45),(D46), 








(D47),(j2l00j), 








(I2.101P 


(D49) = 


[po + + + ed] - So] fd] < p + + + ed] 


-fd] 


by (D15),(D48), 








fl2.34p.fl2.38p 
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(D50) 


Po + + + erf] - 5 ] /a| < io 




by (Premise), (j2.78j) 










(|2. lOOp 


(D51) 


=^ [t - t ] < [x - [[[x + + + erf] - 6 ] 





- fd]] 


by (12.311) 


(D52) 


[x - Xq] < Sj. 






by (C1),(|2.100P 


(D53) 


=>► [1 - 1] < £i 






bv rt2.43p.fl2.10ip 


(Z>54) 


[[1 - 1] + [x + + + ed]] < [d + [t - 


f + + 


erf]] 


by (D2),fl23U) 


(Z>55) 


[[1 - 1] + [x + + + ed]] < [[.t + + - 


1- erf] 4 


*o] 


bv (I2.44p.fl2.100p 


(L>56) 


[[io + + + erf] + [1 - 1]] < [[x + + - 


h erf] 4 


<*o] 


bv fl2.44p.fl2.10ip 


(L>57) 


[[[io + + + erf] + 1] - 1] < [[x + + - 


h erf] 4 


So] 


bv fl2.48p.fl2.10ip 


(£>58) 


=> [[[*o + + + erf] - 1] + 1] < [[£„ + + - 


(- erf] 4 


So] 


bv fl2.46p.fl2.10ip 


(D59) 


[t + + + erf] < [[i + + + erf] 


+ S ] 




bv fl2.45p.fl2.10ip 


(D60) 


=>- [x + + + erf] = 










[po + + + erf] fd] + + + 


erf] 




by (B2),(J2T7D 


(DQ1) 


[[so + + + erf] + 5 ] = 










[p + + + erf] + 5 ] + "+ 


- + erf] 




by (D34),P777I) 


(D62) 


=> [po + + + erf] fd] 4 + + 


erf] < 








[[[[xq 4 4 4 erf] + S ] fd] 4 4 


■4 erf] 




by (D59),(D60), 

(D61),(|2.100P, 

(I2.101P 


(£>63) = 


=» Po + + + erf] /rf] < [p + + + erf] 


4<5 ]- 


fd] 


by (D35),(|2T7|>, 
(D62),(|2.38p 


(D64) 


T < [p + + + erf] + So] 


-fd] 




by (Premise), (J2.78P 
(I2.101P 


(.D65) 


=>- [x - Xq] < \[[[xq 4 + + erf] + Sq] 


-fd] 


— Tnl 


by (12.311) 


(D66) 


[1 - 1] < [[p 4 4 4 erf] 4 S ] 


- fd] - 


UJ 


bv fl2.43P.fl2.10lP 




=>[1-1]<5 2 






bv fC2Ul2.100P 


(D68) 


[1 - 1] < 5 

L J 






bv (D5) (D53) 
(D67) 


(D69) 


=> [[1 - 1] + 1] < [S + 1] 

L L J 1 J L 1 J 






by f]2.3ip 


(L>70) 


=>■ 1< [S 4 1] 






bv (12. 451), (12. 1011) 


(72711 


1< [1 + <S] 






bv (12 441) (12 100P 


(D72) 


1 < T 






by (Premise) 


(D73) 


=Mi + <s] < [*o + 5] 






by (]2.31|) 


(D74) 


^ 1 < [t + 5] 






by (D71),(D73), 
(J2.98P 


(£75) 


=^ [5 + io] < [5 2 + so] 






bv (D7).(|2.31P 
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(D76) 




[x + 5]< [5 2 + x ] 




bv rt2.44p.fl2.10ip 


(D77) 




[x + 8]< [x + S 2 ] 




bv O2.44p.fl2.100p 


(D78) 




[[x + 5] + + + erf) < [[x + 82] 4 


- + + ed] 


by (D74),(D77), 










(]2.37|) 






[[x + 5] + + + ed] < [[z + + + 


■ ed] + e] 


by (D78),(D42), 










(|2.98|) 


(D80) 




=►/([£„ + 5]) < [f(xo) + e 


1 


by (Premise) 






Since /(£) is strictly increasing, f(x — 


5) < /(i) 


and f(x) < f(x + 5) hold 




for all x G (x — 5,x + 5). (D23) and (D80) derive that [f(x ) — e] < f(x) and 




m 


1 < [f( x o) + £ ] hold f° r an ^ e (^0 — 5, x + 5)- 




(2) [[i„ 


+ + + ed] -S ] < 1. 










We construct 5 as follows: 










(£71) 8=[[[[x + + + eJ\ + 8o\- 


fd] 


— Xq] 


(Fl) 




[[x + + + erf] - S ] < 1 




by (Premise) 


(F2) 




[[xo + + + ed] - [e - -[1 + 1]]] < 1 


by (B18),(j2.1()l|) 


(F3) 




[1 - 1] < e 




by (Premise) 


(F4) 




K[l + 1] 




by (j2.30|) 


(F5) 




=► [1 - 1] < [[1 + 1] - 1] 




by (j2.31|) 


(F6) 




=► [1 - 1] < [[1 - 1] + 1] 




by (|2.46|) 


(FT) 




=► [1 - 1] < 1 




by (j2.45|) 


(F8) 




=► [£--[1 + l]]<[e--l] 




by (F3),(F4),(F7),(]2Ji3j) 


(F9) 




[e--l]=e 




by (]2.51|) 


(F10) 




[£--[1 + 1]] < £ 




by (F8),(F9),f|2.100p 


(Fll) 




=> [[io + + + ed] - e] < 










[[x + + + ed]-[e- -[1 + 1]]] 




by (]2.31j) 


(F12) 




[[x + + + ed] - e] < 1 




by (F11),(F2),(|2.98P 


(F13) 




[/(*o) - e] < 1 






(F14) 




/(l) = [1 + + + ed] 






(F15) 




=►/(!) = ! 




by (j2.73|) 


(F16) 




=}► 1 < /(£) 




by (A1)~(A8),(|2.1()1|> 


(F17) 




[/(xo) - e] < /(*) 




by (F13),(F16),(j2.98!) 


(F18) 




«J = [£--[1 + 1]] 




by (B18) 


(F19) 




[1 - 1] < 8 




by (B7),(F18), 










(]2.100P 


(F20) 




[e --[1 + 1]]<£ 




by (D9) 


(F21) 




=^ <5 < £ 




by (F20),(F18), 
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(F22) S=[[[[x + + + ed] + 8 }- -fd]-x ] by (El) 

(F23) [S + ^] = [*o + [[[[^o + + + ed] + <J ] fd] - x )] by Q23HD 

(F24) [x + 5] = [[[[[S + + + ed] + £„] /d] - x ] + S ] by <^M,^HB 

(F25) ^[x + 5] = po + + + ed] + 5 ]- - /d] bv f[2^. (jm3]) 

(F26) [1 - 1] < S by (F19) 

(F27) => [[1 - 1] + 1] < [So + 1] by (F26),d23H> 

(F28) 1 < [cJ + 1] by (|2^5|> . (|2^DIj) 

(F29) 1< [1 + <J ] by (raji. lQDDI) 

(F30) 1 < [x + + + ed] by (B2) 

(F31) [1 + cJ ] < [[So + + + ed] + cJ ] by ([MI]) 

(F32) 1 < [[So + + + ed] + 8 Q ] by (F29),(F31), 

(F33) 1 < [p + + + ed] + 5o] /d] by (TOD 

(F34) =4> [[x + S] + + + ed] = 

[[[[xo + + + ed] + 5 ] fd] + + + ed] by (F25),(F33), 

(F35) ^[[x + S] + + + ed) = [[x Q + + + ed) + 5 ] by (EIZZD , (EZI3D 

(F36) (J < £ by (F21) 

(F37) [So + [So + + + ed]] < [e + [i„ + + + ed]] by (TOT]) 

(F38) [[x + + + ed] + (J ] < [e + [x + + + ed]] by ([2311) , HH 

(F39) [[x + + + ed] + cJ ] < [[S + + + ed] + e] by , flnDDJ 

(F40) [[S + (J] + + + ed] < [[x + + + ed] + £] by (F39),(F35), 

(l^TOl]) 

(F41) =4- /([x + <J]) < [/(So) + e] by (Premise) 

Since f(x) is strictly increasing, /(l) < f(x) and /(x) < f(xo + 5) hold for 
all x G (S — 5, So + S). (F17) and (F41) derive that [/(So) — e] < /(S) and 
/(S) < [/(So) + £ ] hold for all x G (x — S, x + 5). 
Items 2(d)N2(d)ii derive that f(x) is continuous. 

(F16) derives that 1 < f(x) holds on the domain [l,+oo). For any number 1 < e, 

f)2.34p and (I2.74p always derive that [[e fd] + 1] G [1, +oo). So there always exists 

x — [[e fd] + 1] on the domain [1, +oo). 

(Gl) [1 - 1] < 1 by (F7) 

(G2) =*[[l-l] + [e /d]]<[l + [s fd]] hy(gm 

(C73) =>[[e- -/d] + [l-l]]<[l + [e- - fd]] by (g3j,pTI)T|) 

(G4) =►[[[£- -/d] + i]-i]<[i+[ £ - -/d]] by ozasD.azzniiD 

(C75) =►[[[£- -/d]-l] + l]<[l + [e /d]] by i^M),<^m 
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(G6) =>[e fd]<[l + [e fd}} by ([132]), <&M§ 

(G7) =>[e- -/d]<[[£- -M] + l] by (^35),d23DQD 

(G8) x = [[e fd] + 1] by (Premise) 

(G9) =>[e- -/d]<x by (G7),(G8),(J2I00D 

(G10) =*[[e /d] + + + erf] < [x + + + ed] by <^M,^M,<^M 

(Gil) ^£<[x + + + ed) bv fPremise). (l2J7j) . (l2T0B 

(G12) =^ e < f(x ) by (Premise) 



(G1)~(G12) derive that f(x) is unbounded. Since f(x) is both continuous and un- 
bounded on the domain [1, +oo), it maps the domain [1, +oo) to the range [1, +oo). □ 

Theorem 4.4. For a function f : R \ (— oo, [1 — 1]) — >■ R defined by f(x) = [e + + + ex], 
it is injective continuous and maps the domain [[1 — 1], +oo) to the range [1, +oo). 

Proof. According to (I2.19p . the symbol 'e' represents some successive '+' — "+■■•+". 
According to (I2.20p . the symbol '/' represents some successive '— ' — "— According 
to (I2.2ip . the symbol 'g' represents some successive '/' — "/ •••/"■ 



(Al) Suppose that Xi,x 2 G [[1 — l],+oo) with X\ ^ x 2 . 





Without loss of generality, we assume that x\ < x 2 . 




(A2) 


([1 - 1] < ii) A (xi < x 2 ) 




(A3) 


1 < e 


by (Premise) 


(AA) 


=> [e + + + exi] < [e + + + ex 2 } 


by (A2),(A3),(|2.39|) 


(A5) 


f(xi) = [e + + + exi] 


by (Premise) 


(A6) 


=► f(x x ) < [e + + + ex 2 ] 


by (A4),(A5),(]2.101|) 


(A7) 


f(x 2 ) = [e + + + ex 2 ] 


by (Premise) 


(AS) 


=> f{xi) < f(x 2 ) 


by (A6),(A7),(!2100|) 


(A9) 


=► -.(/(ii) = f(x 2 )) 


by (12.97|) 



(A1)~(A9) derive that f(x) is injective. (A1)~(A8) derive that f(x) is strictly increas- 
ing. 

For any number x G [[1 — 1], +oo) and any number [1 — 1] < £, we can always construct 
5 as follows: 



(Bl) [e + + + e[l - 1]] < [e + + + ex ] by $ZM,$ZTM 

(52) =^l<[e + + + ex ] by fl2ZID,GDI2P 

(53) ^ [1 - 1] < [[e + + + ex ] - 1] by (JMU) 

(54) [1 - 1] < e by (Premise) 

(55) =► [[1 - 1] - -[1 + 1]] < [5 - -[1 + 1]] by ^M) 

(B6) =► [[i _ _[i + 1]] _ [i _ _[i + 1]]] <[£--[! + 1]] by ^m, d2nnu) 

(57) [1-1] < [£--[1 + 1]] by (JZ43]), d2ini]) 

(58) =►[[!-!] + [£- -[1 + 1]]] <[[ £ --[l + l]] + [ £ - -[1 + 1]]] by^ 

(B9) [[1- l] + [e- -[1 + 1]]] < [[e + e]--[l + l]] by (J23HD, (J2I00]) 
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(BIO) 


rre — ri h 

LL L 


. lii + n _ iji < rr £ + £ i _ 

JJ 1 L JJ LL 1 J 


L 


+ 111 

1 J J 


bv (I2.44p. (I2.10ip 


(511) 


L L L L 


f 111 + 11 - 11 < \\e + e}- 

1 -LJJ 1 -lj ^LL 1 J 


L 


+ 111 
1 jj 


bv (12.480. (12.1010 


(512) => 


f[[e--[l- 

L LL L 


f 111 - 11 + 11 < \\e + e}- 

i jj ji J L L 1 J 


-fl 


+ 111 

1 J J 


bv (I2.46P. (I2.101P 


(513) 


L 


-[1 + 111 < f[e + el--[l 

L 1 J J L L J L 


+- 111 

1 J J 




bv (12.450. 112.1010 


(514) 


*[e--[l 


+ 1]]<[[ £ + +[1 + 1]]-- 


-[1 + 


111 


bv (I2.55p.fl2.5ip. 












p.ioop 


(515) 


*[e~-[l 


+ l]]<[ e + +[[l + l]-- 


[1 + 


1111 


bv 02.59p.fl2. 100p 


(516) 




[ £ __[1 + 1]] <[ e + + l] 






bv fl2.52p.fl2. 100p 


(517) 




[e- -[1 + 1]] < e 






bv fl2.5ip.fl2.100P 


(518) 




S = [e--[1 + 1]] 








We construct 5 according 


to [[e + + + ex ] - 5 ]. 








(1) l<[[e + 


+ + ex ] - 


5 }. 









We construct 5 as follows: 



(CI) 5i = [x - [[[e + + + ex ] - 5 ]///ge]} 
(C2) 5 2 = [[[[e + + + ex ] + 6 ]///ge] - x ] 
(C3) S = mm{S 1 ,5 2 } 



(51) 






<5o=[e--[l + l]] 


by (B18) 


(52) 






[1 - 1] < S 


by (B7),(D1), 
(I2.100P 


(53) 






[ e __[! + !]] < £ 


by (B17) 


(54) 






<5 < £ 


by (D3),(D1), 
(12.101!) 


(55) 






5 = min{<5i,<5 2 } 


by (C3) 


(56) 






^ < 5i 




(57) 






5 < 5 2 




(58) 




*1 = 


[x -[[[e + + + ex ]-S ]///ge]] 


by (CI) 


(59) - 




-*l] = 


= [x - [x - [[[e + + + ex ] - S ]///ge]]] 


bv (I2.114JI 


(510) = 


[x 


-*l] = 


= [[x - x ] + [[[e + + + ex ] - S ]///ge]] 


bv fl2.50p.fl2.113p 


(511) = 




- *d - 


= [[[[e + + + ex ] - S )///ge\ + [x - x )] 


bv fl2.44p.fl2.113p 


(512) : 


[x 


-Si]- 


= [[[[[e + + + ex ] -S ]///ge] + x ] - x ] 


bv (12.481). 02.1131) 


(513) = 


[x 


-Si]- 


= [[[[[e + + + ex ] -S ]///ge] - x ] + x ] 


bv fl2.46p.fl2.113p 


(514) 






- <Ji] = [[[e + + + ex ] - 5o]///^e] 


bv fl2.45p.(l2.113p 


(515) 




[1- 


1] < [[[e + + + ex ]-6 ]///ge] 


by (Premise), (|2.36| 


(516) 






[1 - 1] < [x - 5i] 


by (D14),(D15), 
(12.1001) 


(517) 






=>- [x - 5i] < [x - 5] 


by (D6).(|2.3ip 
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(D1S) 


=^ fe + + + efx n - < fe + + + efx n - 511 


by (D 17), (|2.39). 






(I2.114P 


(D19) 


=>- [e + + + efffe + + + ex n l - S n )///qe\) < 

' L 1 1 1 LLL U J uj / / / j J J — 






fe + + + efin - 511 


by (D14),(|2l0lj) 


{D20) 


=>- ffe + + + ex l - 5 1 < fe + + + efx - 511 

L |_ 1 u j u j L L J J 


bv (l2.79D.H2.101p 


(D21) 


=>- ffe + + + ex l - e] < ffe + + + ex l - S ) 

LL uj j L L <J j u j 


by IDA), (12.311) 


(D22) 


=>- ffe + + + exnl - el < fe + + + efx - 511 

[_ l_ 1 1 1 uj j L L *-* J J 


by (D21UD20), 






(]2.98| 


{D23) 


\f(x )-e] < f(\x Q -S)) 


by (Premise) 


(D2A) 


£ 2 = [[[[e + + + exnl + <Snl ///gel - x l 

z LLL u j 1 *J J / / / j J uj 


by (C2) 


(D25) 


=>■ fx n + 5ol = fxn + ffffe + + + exnl + 1 Iqe] — Xnll 

L U 1 L U 1 LLLL U J 1 UJ / / / iy J ujj 


bv (12.1141) 


(D26) 


fx n + 5 2 1 = fffffe + + + ex n l + 5o\///qe] - Xn) + Xnl 

' l u 1 £\ LLLLL u j 1 / / iv J ^ u j 1 *^ u j 


bv (12.441). (12. 1131) 


(D27) 


=>• f^o + 5 2 1 = fffe + + + ex n l + Sn]///qe\ 

L U 1 zj LLL U J 1 UJ / / / j J 


bv (12.451). (12. 1131) 


(D28) 


fl - 11 < 5n 
L J u 


by (D2) 


(D29) 


=^ ffl - 11 + 11 < f5n + 11 
LL J J [ u i -*-j 


by (D28),(|2.31|) 


(D30) 


=>- 1 < \Sn + 11 
L u J 


bv (12.451). (12. 1011) 


(D31) 


=>- i < fi + <y 

L u j 


bv (12.441). (12. 1001) 


{D32) 


1 < fe + + + ex n l 

|_— i i i U J 


by (B2) 


(D33) 


=>- fl + 8n\ < ffe + + + einl + So) 

' y -*- i u j LL u j 1 u j 


by (12.311) 


(L>34) 


=>- 1 < ffe + + + exnl + So) 

LL uj 1 uj 


by (D31),(D33), 






(j2.98j) 


(D35) 


=> fl - 11 < fffe + + + exnl + So)///qe) 

i j-j [LL u j 1 uj / / / j j 


by (12.351) 


(D36) 


=>- fe + + + efx n + 5 2 11 = 

L 1 l u 1 z j j 






fe + + + efffe + + + ex l + 5 l///oell 

L 1 1 1 LLL u j 1 uj/ / / iy JJ 


by (D27),(D35), 






(I2.114P 


(D37) 


=>• fe + + + efx n + 5 2 11 = ffe + + + einl + 5 n l 

' i i 1 U 1 J J II ' ' ' — ^ j i u J 


bv (I2.79P.(I2.113P 


(D38) 


5o < B 


by (D4) 


{D39) 


f5 n + fe + + + exnll < \e + fe + + + exn)) 

L u 1 l ujj l L u j j 


by (I2.31P 


(DAO) 


=>• [fe + + + ex ] + 5 ] < fe + [e + + + e£ ]] 


bv (12.441). (12.1011) 


(DAI) 


=>• ffe + + + exnl + 5 1 < ffe + + + ex l + ^1 

LL uj 1 uj LL uj 1 ^j 


bv (12.441). (12. 100p 


(DA2) 


=>• fe + + + efx n + 5 2 11 < ffe + + + ex n l + ^1 

L 1 L " -^JJ LL uj 1 j 


by (D41),(D37), 






(12.1011) 


(D43) 


=^ [fe + + + ex ] - 5 ] < [fe + + + ex ] - [1 - 1]] 


by (D2),(|2Ji]J) 


(L>44) 


[fe + + + ex ] - 5 ] < [ffe + + + ex ] - 1] + 1] 


bv (I2.50I).(I2.100P 


(D45) 


[fe + + + ex ] - 5 ] < [e + + + ef ] 


bv (12.451). (I2.100P 


(DA6) 


=^ [fe + + + ex ] - 5 ] = 






[e + + + efffe + + + ex ] - S )///ge)) 


by (Premise), ((2/79J) 
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(1)47) 
(D48) 



(D4Q) 

(D50) 

(D51) 
(£>52) 
(D53) 
(D54) 
(£>55) 
(£>56) 
(D57) 
(D58) 
(£>59) 
(£>60) 

(£>61) 

(D62) 



(D63) 

(D64) 

(D65) 
(D66) 
(D67) 
(D68) 

(D69) 



=► [e + + + ex ] = 
[e + + + e[[e + + + ex ]///ge]] 
[e + + + e[[[e + + + ex ] - S ]///ge]] < 
[e + + + e[[e + + + ex ]///ge]] 



by (B2),d22i 



by (D45),(D46), 
(D47). (j27TUU|) . 

[e + + + ex ] - 5 ]///ge] < [[e + + + ex ]///ge] by (D15),(D48), 



=► [[[e + + + ex ] - 6 ]///ge] < x 

[xq-Xq] < [x - [[[e + + + ex ] -S ]///ge]] 
=► [x - Xq] < 5i 



< [S + [e + + + ex ]] 

< [[e + + + ex ] + So] 

< [[e + + + ex ] + S ] 

< [[e + + + ex ] + (J ] 

< [[e + + + ex ] + S ] 



=►[1-1 
[[1 - 1] + [e + + + ex ] 
[[1 - 1] + [e + + + ex ] 
[[e + + + ex ] + [1 - 1] 
[[[e + + + ex ] + 1] - 1 
[[[e + + + ex ] - 1] + 1 
=► [e + + + ex ] < [[e + + + ex ] + 5 
=► [e + + + ex ] = 
[e + + + e[[e + + + ex ]///^e]] 
=► [[e + + + exo] + So] = 
[e + + + e[[[e + + + ex ] + S ]///ge]] 
=► [e + + + e[[e + + + ex ]///ge]] < 
[e + + + e[[[e + + + ex ] + 6o]///ge\] 



(JI35D,(J230D 

by (Premise), ([220]), 

(MID 

by (E3I]) 

by (Cl), d23DQP 
by d233p, p3n0 

by (D2),H23U 
by (127441). (1271001) 
by (12714]) . (127107]) 

by (CT. CTO 
by (12746]) . (12~T0T]) 
by (12743]) . (127T0T]) 

by (B2),fl2Zi 

by (D34),d22SD 



by (D59),(D60), 

(D6i),d2nnnD, 

(127TDT]) 

[[e + + + ex ]///ge] < [[[e + + + ex ] + <J ]///^e] by (D35),(D62), 



=► x < [[[e + + + ex ] + 6 ]///ge] 

[xo - x ] < [[[[e + + + ex ] + S ]///ge] - x ] 
> [1 - 1] < [[[[e + + + ex ] + 5 ]///ge] - x ] 
=► [1 - 1] < (J 2 
[1 - 1] < 5 

[1 - 1] < x 



by (Premise), Q2ISQJ), 

fl^TDT]) 

by (E3U) 

by (12743]) . (12~T0T]) 

by (C2),o2nnn]) 

by (D5),(D53), 
(D67) 

by (Premise) 
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(D70) 


=> [fl -1] + S\< [i + 5] 

LL J ' J — L U ' J 


by (|2.31|) 


(D71) 


=>- [5 + fl - 111 < fen + 51 

L L J J L ^ J 


bv (l2.44D.H2.10ip 


(D72) 


[[5 + 11 - 11 < ten + 51 

LL J J — L VJ 1 J 


bv (12.481). (I2.101P 


(D73) 


\\S - 11 + 11 < fin + 51 

' LL J J — L u 1 


bv (12.461). (12. 1011) 


(D7A) 


=>- 5 < ten + 51 

L *J J 


bv (12.451). (12. 10ip 


(D75) 


=>» fl - 11 < ten + 51 

L J L " J 


by (D68),(D74), 






(12.981) 

M — : y 


(D76) 


=>- [5 + Xol < \5o + Xn] 


by (D7), (12.311) 


(D77) 


=>- ten + 51 < [5 2 + inl 

L " J — L ** u j 


bv (12.441). (12.1011) 

j \j 1/ 1 \) y 


(D78) 


=>• ten + 51 < ten + 69) 


bv (12 441) (12 1001) 


(D79) 


[e + + + e[x + 5]] < [e + + + e[x + 5 2 ]] 


by (D75),(D78), 






(12.391) 


(D80) 


[e + + + e[x + 5]] < [fe + + + ei ] + e] 


by (D79),(D42), 






02.980 


(D81) 


/([*„ + <*]) < [f(x ) + e] 


by (Premise) 



Since /(x) is strictly increasing, /(xo — 5) < /(i) and /(i) < /(io + 5) hold 
for all x G (io — 5, x + 5). (D23) and (D81) derive that [/(io) — £ ] < /(i) an d 
/(i) < [/(io) + e] hold for all i G (io — 5, x + 5). 
(2) [[e + + + ex ] - 5 ] < 1. 



We construct 5 as follows: 

(Fl) 5 = [[[[e + + + ex ] + 6 ]///ge] - i ] 



(Fl) 


[[e + + + ex ] - S ] < 1 


by (Premise) 


(F2) 


=> [[e + + + ex ] - [e - -[1 + 1]]] < 1 


by (B18),(]2.101|) 


(F3) 


[1 - 1] < 8 


by (Premise) 


(F4) 


K [1 + 1] 


by (|2.30| 


(F5) 


=► [1 - 1] < [[1 + 1] - 1] 


by (j2.31|) 


(F6) 


=► [1 - 1] < [[1 - 1] + 1] 


by (|2.46|) 


(F7) 


=► [1 - 1] < 1 


bv (I2.45|) 


(F8) 


=► [e--[l + l]]<[e--l] 


by (F3), (14), (17), ((2,331) 


(F9) 


[e--l]=e 


by (j2.51|) 


(F10) 


[£--[1 + 1]] <e 


by (F8),(F9),(]2.100P 


(Fll) 


[[e + + + ex ] -e] < 






[[e + + + ex ] - [e- -[1 + 1]]] 


bv (l2.3l|) 


(F12) 


[[e + + + ei ] - e] < 1 


by (111),(F2),(]2.98P 


(F13) 


[/(io) - e] < 1 




(Fll) 


/([l-l]) = [e + + + e[l-l]] 




(F15) 


=►/([! -1]) = 1 


by (j2.71|) 
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(F16) =► 1 < f(x) by (Al)~(A8),Q2inU 

(F17) [/(*„) - e] < /(£) by (F13),(F16),(JI98D 

(F18) So = [e--[l + l]\ by(B18) 

(F19) [1 - 1] < tfo by (B7),(F18), 

( EDQD 

(F20) [e--[l + l]]<e by (D9) 

(F21) =>• <J < e by (F20),(F18), 

( 1230T]) 

(F22) 5 = [[[[e + + + ex ] + S Q ]///ge] - x ] by (El) 

(F23) ^[xo + 8] = [x + [[[[e + + + ex }+6o]///ge}-xo\] bv f l2TT4|) 

(F24) ^[x + (J] = [[[[[e + + + exo] + (Jo]///^]-x ]+Xo] by 023^,(123135 

(F25) ^[x + 5] = [[[e + + + ex ] + 5 ]///ge] by (JZ35J) , (J23I3]) 

(F26) [1 - 1] < <J by (F19) 

(F27) [[1 - 1] + 1] < [(J + 1] by (F26),d23B 

(F28) 1 < [(J + 1] by (J23SJ) , (EHHID 

(F29) 1< [1 + (Jo] by flCTj).f l2TT00|) 

(F30) 1 < [e + + + ex ] by (B2) 

(F31) [1 + (J ] < [[e + + + ex ] + <J ] by fl23D 

(F32) 1< [[e + + + ex ] + <J ] by (F29),(F31), 

(F33) [1-1] < [[[e + + + ex ]+S ]///ge] by <^M> 

(F34) [e + + + e[x + 5]] = 

[e + + + e[[[e + + + ex ] + S ]///ge]} by (F25),(F33), 

( EED 

(F35) [e + + + e[x + 5]] = [[e + + + ex ] + S ] by (J2ZZSI> , dHHBH 

(F36) (J < £ by (F21) 

(F37) [So + [e + + + ex ]] < [e + [e + + + ex ]] by Q23D 

(F38) [[e + + + ex ] + 5 ] < [e + [e + + + ex ]] by (1231, (EM) 

(F39) [[e + + + ex ] + <J ] < [[e + + + ex ] + e] by (J23HI , (I23DQD 

(F40) => [e + + + e[x + <J]] < [[e + + + ex ] +e] by (F39),(F35), 

(12301]) 

(F41) /([x + <J]) < [/(x ) + e] by (Premise) 



Since /(x) is strictly increasing, /(l) < /(x) and /(x) < /(xq + S) hold for 
all x £ (xo — 5, io + 5). (F17) and (F41) derive that [/(xo) — e] < /(x) and 
/(x) < [/(xq) + e] hold for all x £ (x — S,Xq + S). 



Items 2(d)N2(d)ii derive that /(x) is continuous. 
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(F16) derives that 1 < f(x) holds on the domain [l,+oo). For any number 1 < e, 
(I2.35P and (I2.75P always derive that [[e/ //ge] + 1] G [[1 — 1], +oo). So there always exists 
xo = [[e/ / /ge] + 1] on the domain [1, +oo). 

(Gl) [1 - 1] < 1 by (F7) 

(G2) ^[[l-l] + [s///ge]]<[l + [s///ge\] by (£3U) 

(G3) ^[[e///ge] + [l-l]]<[l + [e///ge]] by 

(G4) =*► [[[e///</e] + 1] - 1] < [1 + [e///<7e]] by (gaBXplOl]) 

(G5) ^[[[e///ge}-l] + l]<[l + [e///ge]] by (fOE) ,02201) 

(G6) [e/J/ge] < [1 + [e///^]] by (EZ5D , jM} 

(G7) [e///(/e] < [[e///^e] + 1] by ggg^piPD) 

(G8) x = [[e///ge] + 1] by (Premise) 

(G9) [ £ ///<?e] < x by (G7),(G8),@3DD) 

(G10) [e + + + e[e///ge]} < [e + + + ex ] by , , (EMD 

(Gil) =^£< [e + + + ex ] bv fPremise). f[2Tm f l23011 

(G12) =>• e < f(x ) by (Premise) 

(G1)~(G12) derive that f(x) is unbounded. Since f(x) is both continuous and un- 
bounded on the domain [[1 — l],+oo), it maps the domain [[1 — l],+oo) to the range 
[l,+oo). □ 

In the following, we classify the binary numerical computations for high operations. 

(1) Binary numerical computation for operations + + ++, , ////. 

(a) A binary numerical computation for [[ax a 2 ] + + + +n] with 1 < [&i a 2 ] 

can be achieved. 

{AI) [[at - -a 2 ] + + + +n] = 

[[&! - -a 2 ] + + + [[ax - -a 2 ] + + + +[h - 1]]] by f[2~8Tj) 

(A2) Let us distinguish these [ax a 2 ] with 

the subscripts {(1), (2), (3), • • • } 
(A3) [[ai - -a 2 ] + + + +n] = 

[[ai a 2 ] ( i) + + + [[ai a 2 ] {2) + + H 

[[ai - -a 2 ](n-i) + + + [[ai ~ -«2](n) + + + +1]] • • ■]] by 
(A4) [[ai - -a 2 ] + + + +n] = 

[[ai a 2 ] ( i) + + + [[ai a 2 ] (2) + + H 

[[ai - -a 2 ](n-i) + + + [ai- -a 2 ] ( n)] • • •]] by ([220]) 



Then (Al)~(A4) have reduced one + + ++ operation to many + + + 
operations. Since §4.41 has achieved the binary numerical computation for 

[[fii a 2 ] + + + 72.], (Al)~(A4) can achieve a binary numerical computation 

for [[ax a 2 ] + + + +n]. 
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(b) A binary numerical computation for [[di d 2 ] ft] with 1 < [di d 2 ] 

can be achieved. 

The binary numerical computation for [[di d 2 ] n] is equated 

with the numerical root finding of the equation x = [[d± d 2 ] n]- 

(Al) x = [[ai d 2 ] ft] 



[[[ai - -d 2 ] ft] + + + +n] by (gUj]) 

(A3) [x + + + +n] = [ai - -d 2 ] by (J2T7D , $ZM) , dffl) 

According to fl2.37pMI2.40p . the function f(x) = [x + + + +n] is defined 
on the domain [l,+oo). Theorem 14.31 implies that f(x) is continuous on the 

domain [1, [ai a 2 ]]. 

f l2~T3j) derives that [l++++n] = 1. So f!2.10ip derives that [l++++n] < 

[di a 2 ]. f)2.70p derives that [[di a 2 ] + + + + l] = [di a 2 ]. f )2.39p derives 

that [[di a 2 ] + + + +1] < [[a x a 2 ] + + + +n}. So (I2.1l)ip derives that 

[di d 2 ] < [[&i a 2 ] + + + +n]. In summary, both [l + + + +n] < [di a 2 ] 

and [d\ a 2 ] < [[di a 2 ] + + + +n] hold. 

Then Intermediate Value Theorem derives that f(x) has only one root 

on the domain [1, [d\ a 2 ]]. Since Theorem 14.31 implies that f(x) has no 

root on the domain ([d\ a 2 ], +00), f(x) has only one root on the domain 

[l,+oo). Since f(x) belongs to the special equation in §31 Brent's method 
can find the only root of f(x) and achieve the binary numerical computation 

for [[dx d 2 ] ft]- 

(c) A binary numerical computation for [[di a 2 ] + + + +[bi 6 2 ]] with 

1 < [di a 2 ] can be achieved. 



(A2) 



=^ [x + + + +n] 



(i) [b 



b 2 ] < 1. 



(AI) 



I([b 
[T_ci_c 2 . 



M) = 

-±_Ci_C 2 . 



by (EA), 



(A2) 



=> [[Si - -d 2 ] + + + +[b 1 - -b 2 ]\ 
[[ai - -d 2 ] + + + + 



(A3) 



[T_ci_c 2 _ -L_ci_c 2 _]] 

=> [[Si - -02] + + + +[&!- -& 2 ]] 
[[[di d 2 ] + + + +T_ci_c 2 _] 



by (EHU) 



±_Ci_C 2 . 



by (l2T90|).f l2TT3|) 



The variable [T_ci_c 2 _ -L_Ci_c 2 _] represents some irreducible 

fraction in the Farey sequence. 61,62 S N can derive that T_ci_c 2 _, 
±_ci_c 2 _ G A r . Then the items [la] and [lb] can further achieve the 
binary numerical computation for (A3). 
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(ii) 1 < [&i b 2 ). 

(51) [[k - -b 2 }\ = n 

(52) [[h - -b 2 ] - n] < 1 

(53) I([[b 1 --b 2 ]-n}) = 

[T_C!_c 2 _- -±_C!_c 2 _] by (B1),(B2), 

(EA), 

(54) =► [[at - -S 2 ] + + + +[k - -k)} = 

[[fii S 2 ] + + + [[Si S 2 ] 

+ + + + [[6i - -h] - 1]]] by dMU) 

(55) Let us distinguish these [Si S 2 ] 

with the subscripts 
{(1),(2),(3),..-} 

(56) =>> [[ai- -S 2 ] + + + +[&!- -b 2 }} = 

[[Si a 2 ](i) + + + [[Si a 2 ]( 2 ) 

+ + H [[Si a 2 ](n-i) + ++ 

[[Ol - -02](n) + + + + [[&1 - -6 3 ] - 

••■]] by ([MID 

(57) [[Si - -S 2 ] + + + +[6i - -63]] = 
[[Si a 2 ] ( i) + + + [[&i a 2 ](2) 

+ + H [[Si S 2 ]( ft _i) + ++ 

[[Si fi 2 ](ft) + + + + 

[T_C!_c 2 _ - -i-_ci_c 2 _]]] ■ ■ ■ ]] by (B3),dZH^ 



Then the item 2(c)i can further achieve the binary numerical com- 
putation for (B7). 

(d) A binary numerical computation for [[Si S 2 ] [b\ b 2 ]] with 

1 < [Si S 2 ] can be achieved. 

The binary numerical computation for [[Si S 2 ] [61 6 2 ]] is 

equated with the numerical root finding of the equation x = [[Si S 2 ] — 

[&! 6 2 ]]. 

{AI) i=[[Si--S 2 ] [61 --6 2 ]] 

(A2) =4- [x + + + +[h - -b 2 ]] = 

[[[Si S 2 ] [61 62]] 

+ + + + [61 - -62]] by fl2TTl4j) 

(A3) =^[x + + + +[b 1 --b 2 ]] = [ai--<k] by (jT77]), 0221, 

(i) 1< [61 6 3 ]- 
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According to ([Q7) ~ (j2l7j]) , the function f(x) = [x + + + + [h - 
—62]] is defined on the domain [1, +00). Theorem 14.31 implies that f(x) 
is continuous on the domain [1, [hi d 2 ]]. 

( 1X731) derives that [1 + + + +[b x - -b 2 ]} = 1. So (1XTUT]) derives 

that [1 + + ++[foi b 2 ]] < [ax d 2 \. (EZZOD derives that [[fii a 2 ] + 

+ + +1] = [a x fi 2 ]. Since 1 < [fii fi 2 ], (ESHD and ()2Z3])derive 

that [[ax fi 2 ] + + + +1] < [[fi x fi 2 ] + + + +[h b 2 ]\. So 

f!2.10ip derives that [ax d 2 ] < [[ax fi 2 ] + + + +[bx 62]]- In 

summary, both [1 + + + +[&i b 2 ]] < [ax d 2 ] and [dx a 2 ] < 

[[ax d 2 ] + + + +[bx b 2 \] hold. 

Then Intermediate Value Theorem derives that f(x) has only one 

root on the domain [1, [dx 02]]- Since Theorem 14.31 implies that f(x) 

has no root on the domain ([dx d 2 ], +00), f(x) has only one root on 

the domain [l,+oo). Since f(x) belongs to the special equation in §31 
Brent's method can find the only root of f(x) and achieve the binary 

numerical computation for [[dx d 2 ] [61 b 2 }}. 

(ii) [bx - -b 2 \ < 1. 

(51) I{[bx - -b 2 \) = 

[T_ci_c 2 _ - -±_ci_c 2 _] by (EA), 

(ESSD-dUl 

(52) [x + + + -62]] = 

[x + + + + [T_ci_c 2 _ - - ±_cx_c 2 _]] by ( l2~TT4j) 

(53) [x + + + +[bx - -b 2 }} = 

[[x + + + +t_c!_c 2 j ±_ci_c 2 j by ffzm. fium 

(54) =>■ [[x + + + +T_ci_c 2 _] ^_ci_c 2 _] = 

[dx - -d 2 ] by (A3), ([CT3D 

(55) [[fii - -fi 2 ] + + + +[1 - 1]] < 

[[dx - -d 2 ] + + + +±_cx_c 2 _] by flQ9|l 

(B6) => 1 < [[dx - -d 2 ] + + + +±_C!_c 2 _] by fl2TT]).f l23uT]) 

(Bl) ([[dx - -a 2 ] + + + +±_ci_c 2 _]) 

(B8) [[[x + + + +T_cx_c 2 _] ±_ci_c 2 _] 

+ + + + ^_ci_c 2 _] = [[fii fi 2 ] + + + + 

^_ci_c 2 _] by (B4),(B7), 

(59) [x + + + +T_ci_c 2 _] = [[Si fi 2 ] + + + + 

^_ci_c 2 _] by (12771) .( 12TT3D 

(510) 1 < [[[dx - -fi 2 ] + + + +±_C!_c 2 _] 

T_C!_c 2 _] by (B6),d231, 
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(511) ([[[hi - -a 2 ] + + + +±_c 1 _c 2 _] by (H99]) 

T_ci_c 2 _]) 

(512) [[x + + + +T_C!_c 2 _] T_C!_c 2 _] = 

[[[ai a 2 ] + + + +-L_ci_c 2 _] 

T_cx_c 2 _] by (B9),(B11), 

(EHUD 

(513) x = [[[ai h 2 ] + + + +±_ci_c 2 _] 

T_ci_c 2 _] by lEZED , (AI 14) 

The variable [T_Ci_c 2 _ -L_ci_c 2 _] represents some irreducible 

fraction in the Farey sequence. 61,62 6 N can derive that T_ci_c 2 _, 
_L_ci_c 2 _ G iV. Then the items [la] and [lb] can further achieve the 
binary numerical computation for (B13). 

(e) A binary numerical computation for [[hi h 2 }/ / / /[bi b 2 ]} with 1 < 

[61 h 2 ] and with 1 < [61 b 2 ] can be achieved. 

The binary numerical computation for [[hi a 2 }/ / / /[bi 62]] is equated 

with the numerical root finding of the equation x = [[hi h 2 ]/ / / /[bi 62 1 



(AI) x = [[hi - -h 2 ]////[bi - -b 2 \] 

(A2) [1 - 1] < [[hi - -h 2 ]////[bi - -b 2 ]] by ([M 

(A3) 1 < [[61 - -6 3 ] + + + + 

[[ai - -a 2 }////[h - -6 2 ]]] by ^M,^ 

(AA) = ^([[& 1 __fc 2 ] + + + + 

[[fii " -a 2 ]////[bi - -6a]]]) by (JZSSJ 

(A5) [[61 - -63] + + + +x] = 

[[h - -6 2 ] + + + + 

- -« 2 ]////[&i - -62]]] by mm 

(A6) [[bi - -6 2 ] + + + +x] = [hi - -h 2 ] by (J2T9]),(|2TI3D 

According to (EE37|l ~( j21fip . the function /(x) = [[61 6 2 ] + + + +x] is 

defined on the domain [[1 — 1], +00). Theorem 14.41 implies that we can always 

increase xq until [hi h 2 ] < [[bi 62] + + + +xo] holds. Theorem 14.41 

implies that f(x) is continuous on the domain [[1 — l],xo]. 

fl27TTj) derives that [[61 b 2 ] + + + +[1 - 1]] = 1. So fl^TOT]) derives 

that [[61 b 2 ] + + + +[1 — 1]] < [hi h 2 }. In summary, both [[61 b 2 ] + 

+ + + [1 - 1]] < [ai h 2 ] and [h t a 2 ] < [[hi b 2 ] + + + +x ] hold. 

Then Intermediate Value Theorem derives that f(x) has only one root 
on the domain [[1 — l],Xo]. Since Theorem 14.41 implies that f(x) has no root 
on the domain (x , +00), f(x) has only one root on the domain [[1 — 1], +00). 
Since f(x) belongs to the special equation in §31 Brent's method can find the 
only root of f(x) and achieve the binary numerical computation for [[hi — 
-h 2 \////[bi--b 2 \]. 
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(2) If the binary numerical computations for [d + + + eb], [a fb], [a// / gb] are 

achieved, then the binary numerical computations for [d + + + +e6], [d fb], 

[a/ 1 1 1 gb] are also achieved. 

According to ( I2.19p . the symbol 'e' represents some successive '+' — "+ •••+". 
According to (I2.20j) . the symbol '/' represents some successive ' — ' — 
According to (12.211) . the symbol 'g J represents some successive '/' — "/ •••/"• 

(a) A binary numerical computation for [[di d 2 ] + + + +efi] with 1 < [Si S 2 ] 

can be achieved. 

(Al) [[di - -d 2 ] + + + +en] = 

[[Si d 2 ] + + + e[[di d 2 ] + + + +e[n - 1]]] by OI) 

(A2) Let us distinguish these [di d 2 ] with 

the subscripts {(1), (2), (3), • • • } 

(A3) =>- [[di - -d 2 ] + + + +en] = 

[[di d 2 ] ( i) + + + e[[di d 2 ] (2 ) + + + e • ■ • 

[[Si - -a 2 ](n-i) + + + e[[fii - -d 2 ] w + + + +el]] • • •]] by (JSP 

(A4) [[di - -d 2 ] + + + +en] = 

[[di d 2 ] ( i) + + + e[[di d 2 ]( 2 ) + + + €■■■ 

[[di d 2 ](n-i) + + + e[di a 2 ](n)] • • ■]] by (|2.70|) 

Then (A1)~(A4) have reduced one + + + + e operation to many + + +e 

operations. Since the binary numerical computation for [[di d 2 ] + + + 

en] is supposed to be achieved, (Al)~(A4) can achieve a binary numerical 
computation for [[di d 2 ] + + + +en]. 

(b) A binary numerical computation for [[di d 2 ] fn] with 1 < [di d 2 ] 

can be achieved. 

The binary numerical computation for [[di d 2 ] fn] is equated 

with the numerical root finding of the equation x = [[di d 2 ] fn]. 

(Al) x = [[di d 2 ] fn] 

(A2) [x + + + +en] = 

[[[di - -d 2 ] fn] + + + +eh] by (I2~TT4|1 

(A3) => [x + + + +en] = [d x - -d 2 ] by ([2177]), ([221, (O 

According to (I2.37p ~( f2.40p . the function f(x) = [x + + + +eh] is defined 
on the domain [l,+oo). Theorem 14.31 implies that f(x) is continuous on the 
domain [1, [di d 2 ]]. 

( 1277511 derives that [1 + + + +eh] = 1. So ( l2.1Uip derives that [1 + + + 

+en] < [di d 2 ]. (I2.70p derives that [[di d 2 ] + + + +el] = [di d 2 ]. 

(12739]) derives that [[di d 2 ] + + + +el] < [[Si S 2 ] + + + +en]. So 

(1 2 . 1 1 j) derives that [di d 2 ] < [[di d 2 ] + + + +eh]. In summary, both 

[1 + + + +eh] < [di d 2 ] and [di d 2 ] < [[di d 2 ] + + + +en] hold. 

Then Intermediate Value Theorem derives that f(x) has only one root 
on the domain [1, [di d 2 ]]. Since Theorem 14.31 implies that f(x) has no 
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root on the domain ([dx a 2 ], +00), f(x) has only one root on the domain 

[l,+oo). Since f(x) belongs to the special equation in §3] Brent's method 
can find the only root of f(x) and achieve the binary numerical computation 
for [[hi d 2 ] fn}. 

A binary numerical computation for [[dx d 2 ] + + + +e[bx b 2 ]] with 

1 < [cii d 2 ] can be achieved. 

(i) [h - -b 2 ] < 1. 

(AI) - -fia]) = 

[T_C!_c 2 _ - -±_C!_c 2 _] by (EA), 

(]2~82l) —02789]) 

(A2) [[ai - -d 2 ] + + + +e[6i - -62]] = 

[[61 a 2 ] + + + +e 

[T_C!_c 2 _ - -±_C!_c 2 J] by (JZUID 

(A3) [[ai - -a 2 ] + + + +e[6i - -6 2 ]] = 

[[[ai a 2 ] + + + +eT_ci_c 2 _] 

/±_C!_c 2 _] by (ESOJ) , (J2H3D 

The variable [T_ci_c 2 _ -L_ci_c 2 _] represents some irreducible 

fraction in the Farey sequence. bx,b 2 E N can derive that T_ci_c 2 _, 
_L_ci_c 2 _ G iV. Then the items [2a] and [2b] can further achieve the 
binary numerical computation for (A3). 

(ii) 1< [h - -b 2 \. 

(Bl) [[k - -b 2 }\ = n 

(B2) [[k - -b 2 ] - h] < 1 

(B3) ^I([[b 1 --b 2 ]-h]) = 

[T_C!_c 2 _ - -±_C!_c 2 _] by (B1),(B2), 

(EA), 

(B4) [[ax d 2 ] + + + +e[bi 62]] = 

[[ax d 2 ] + + + e[[dx d 2 ] 

+ + + + e[[bx - -b 2 ] - 1]]] by (ESI]) 

(B5) Let us distinguish these [ax d 2 ] 

with the subscripts 
{(1),(2),(3),--.} 
(B6) [[dx - -d 2 ] + + + +e[bx - -b 2 }} = 

[[dx a 2 ](i) + + + e[[dx a 2 ] (2) 

+ + +e • • • [[ax a 2 ](n_i) + + + e 

[[dx - -a 2 ]( ft) + + + +e[[6i - -6 2 ] - n]]] 

•■■]] by d2SU> 
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(B7) [[fii - -d 2 ] + + + +e[6i - -62]] = 

[[&! a 2 ](i) + + + e[[ai a 2 ]( 2 ) 

+ + +e ■ • • [[ai a 2 ](n-i) + + + e 

[[ai a 2 ]( ft ) + + + +e 

[T_C!_c 2 _ - -±_C!_c 2 _]]] ■ ■ ■ ]] by (B3),@HD 



Then the item 2(c)i can further achieve the binary numerical com- 
putation for (B7). 

(d) A binary numerical computation for [[di d 2 ] f[b\ 6 2 ]] with 

1 < [ai d 2 ] can be achieved. 

The binary numerical computation for [[d\ d 2 ] 6 2 ]] 

is equated with the numerical root finding of the equation x = [[di d 2 ] — 

f[k--b 2 }}. 

(AI) x = [[&! - -d 2 ] /[6i - -b 2 ]] 

(A2) [x + + + +e[6 1 --62]] = 

[[[fii - -d 2 ] /[&! - -& 2 ]] 

+ + + + e[6i - -6 2 ]] by (gJHD 

(A3) [x + + + +e[6i - -6 2 ]] = [ai - -a 2 ] by (J2Z7]) , (JOID , d22S]) 

(i) 1 < [61 6 2 ] ■ 

According to f[2T3Tj) ~ f l2T40|) . the function /(£) = [x + + + +e[bi - 
—62]] is defined on the domain [1, +00). Theorem 14.31 implies that f(x) 
is continuous on the domain [1, [d\ a 2 ]]. 

(12773] derives that [1 + + + +e[6i - -6 2 ]] = 1. So fl27TQl]) derives 

that [l + + + +e[b 1 b 2 }} < [&i a 2 ]. (LTTOj) derives that [[ai a 2 ] + 

+ + + e i] = [a x a 2 ]. Since 1 < [a a a 2 ], fl2739|) and (ETTpderive 

that [[ai a 2 ] + + + +el] < [[ai a 2 ] + + + +e[6i b 2 }}. So 

f)2.10ip derives that [&i d 2 ] < [[«i d 2 ] + + + +e[6i b 2 \\. In 

summary, both [1 + + + +e[b\ b 2 ]] < [dx d 2 ] and [di d 2 ] < 

[[ai d 2 ] + + + +e[bi b 2 \] hold. 

Then Intermediate Value Theorem derives that f(x) has only one 

root on the domain [1, [di d 2 ]]. Since Theorem 14.31 implies that f(x) 

has no root on the domain ([di d 2 ], +00), f(x) has only one root on 

the domain [1, +00). Since f(x) belongs to the special equation in ^J3j 
Brent's method can find the only root of f(x) and achieve the binary 
numerical computation for [[di d 2 ] f[b\ b 2 )]. 

(ii) [bi - -b 2 ] < 1. 

(Bl) I([b t - -6 2 ]) = 

[T_ci_c 2 _ - -±_ci_c 2 _] by (EA), 

(j2T82])~(r2789]) 

(B2) [x + + + +e[b 1 - -b 2 }] = 

[x + + + +e[T_c 1 _c 2 _ - -±_c 1 _c 2 J\ by fiHTQ 
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(53) [x + + + +e[h - -b 2 }\ = 

[[x + + + +eT_C!_c 2 _] A_ci_c 2 _] by fl^90|) . f l2TT3|) 

(54) [[£ + + + +eT_ci_c 2 _] A_ci_c 2 _] = 

[&! - -d 2 ] by (A3). (im^l 

(55) [[ai - -d 2 ] + + + +e[l - 1]] < 

[[ai - -a 2 ] + + + +e±_ci_c 2 _] by fl^BTij) 

(56) 1 < [[fii - -d 2 ] + + + +el__C!_c 2 _] by (|2TT]l. (l230T]l 

(57) => ([[Si - -d 2 ] + + + +e±_ci_c 2 _]) 

(58) [[[x + + + +eT_ci_c 2 _] A_ci_c 2 _] 

+ + + + e_L_ci_c 2 _] = [[di d 2 ] + + + +e 

^i_c 2 _] by (B4),(B7), 

AMU 

(59) =>• [x + + + +eT_ci_c 2 _] = [[di - -d 2 ] + + + +e 

^_ci_c 2 _] by (gT7D,( EISD 

(510) 1 < [[[&! - -d 2 ] + + + +e±_ci_c 2 _] 

/T_c a _c 2 J by (B6),(H33|I, 

(jZTJD 

(511) ^([[[d 1 --d 2 ] + + + +ei__c 1 _c 2 _] bv flCTj) 

/T_C!_C 2 _]) 

(512) [[x + + + +eT_ci_c 2 _] /T_ci_c 2 _] = 

[[[di d 2 ] + + + +e±_ci_c 2 _] 

/T_C!_c 2 _] by (B9),(B11), 

(513) x = [[[di d 2 ] + + + +e±_C!_c 2 _] 

/T_C!_c 2 _] by d23SH,(AH4) 

The variable [T_ci_c 2 _ -L_ci_c 2 _] represents some irreducible 

fraction in the Farey sequence. bi,b 2 E N can derive that T_ci_c 2 _, 
_L_ci_c 2 _ G iV. Then the items [2a] and [2b] can further achieve the 
binary numerical computation for (B13). 

(e) A binary numerical computation for [[di d 2 ]/ '/ / '/ g[bi b 2 ]} with 1 < 

[di d 2 ] and with 1 < [&i b 2 ] can be achieved. 

The binary numerical computation for [[di a 2 ]/ / / / g[bi b 2 ]] is 

equated with the numerical root finding of the equation x = [[di d 2 ]/ / / /g[br 

-HI 

(Al) x = [[d! - -a.y/f/glk - -b 2 }} 

(A2) [1 - 1] < [[dx - -a 2 ]////y[6i - -b 2 }] by (^M) 

(A3) 1 < [[k - -b 2 ] + + + +e 

[[Si - -Zh]////g[i>i - -62]]] by d2^SD,d2ZZTI> 
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(AA) ([[&! --6 2 ] + + + +e 

- -&2]////9[k - -h]]]) by (E39]) 

(A5) [[&i - -b 2 ] + + + +ex] = 
[[61 62] + + + +e 

[[fix - -&»]// //gtin - -h]}} by jSHJ 

(A6) =>[[h 1 --h 2 ] + + + +ex] = [a 1 --a 2 ] by (J2T9]) , (|2TI3D 

According to ( ET57|) ~ (I2730|) . the function /(£) = [[61 6 2 ] + + + +ex ] is 

defined on the domain [[1 — 1], +00). Theorem 14.41 implies that we can always 

increase x until [hi a 2 ] < [[61 b 2 ] + + + +ex ] holds. Theorem 14.41 

implies that f(x) is continuous on the domain [[1 — l],£o]- 

( EpTj derives that [[h - -b 2 ] + + + +e[l - 1]] = 1. So ( I2JTJT]) derives 

that [[61 b 2 ] + + + +e[l — 1]] < [a x a 2 ]. In summary, both [[bi b 2 ] + 

+ + +e[I - I]] < [hi h 2 ] and [hi h 2 ] < [[hi b 2 ] + + + +ex ] hold. 

Then Intermediate Value Theorem derives that f(x) has only one root 
on the domain [[1 — l],Xo]- Since Theorem 14.41 implies that f(x) has no root 
on the domain (x , +00), f(x) has only one root on the domain [[1 — 1], +00). 
Since f(x) belongs to the special equation in §31 Brent's method can find the 
only root of f(x) and achieve the binary numerical computation for [[hi — 
-h 2 ]////g[bi--b 2 ]]. 

(3) By induction, the binary numerical computations for [a + + + +b], [d + + + + + b], 

[a + + + + + +b], [a 6], [a 6], [a 6], [h////b], 

[hj J / / /b], [a///// fb], ■ ■ ■ are all achieved. 
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